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We construct a model for the quark-gluon vertex of Landau gauge QCD. This is of twofold interest: on the one 
hand the quark-gluon interaction is at the heart of quark confinement, on the other hand it is a central element 
in hadron phenomenology based on QCD Greens functions. We employ the non-Abelian one-loop diagram in 
perturbation theory, which is of order N a . As a novelty we replace the tree- level quark and gluon propagators in 
this diagram by their dressed counterparts solving the Dyson-Schwinger equations. The iV c -suppressed Abelian 
diagram is an order of magnitude smaller in various kinematics. We also study the effect of ghost dressing factors 
on the vertex obtaining a construction in good agreement with recent low-momentum lattice calculations. 



1. Construction of the vertex model 

The infrared suppression of the gluon 2-point 
function in QCD |1I2| entails that the bare qqg 
vertex, usually employed in the rainbow trun- 
cation of the Dyson-Schwinger Equation (DSE), 
is insufficient to trigger dynamical chiral sym- 
metry breaking. Therefore we expect an in- 
frared enhancement in the quark and gluon ver- 
tex, as suggested also by its Slavnov- Taylor Iden- 
tity (STI). A model vertex with a Ball-Chiu 
or Curtis-Pennington structure multiplied by en- 
hancing ghost factors has been successfully em- 
ployed in [3] in the quark DSE. For a short sum- 
mary of these results see ref. @]. 

In this work we report a diagrammatic con- 
struction based on the one-loop perturbative 
QCD corrections to the bare vertex 

There are two relevant diagrams, to which we 
refer as Abelian and non-Abelian due to the ver- 
tex attached to the gluon (as is customary). In 
both diagrams we substitute the quark and gluon 
propagators by their dressed counterparts solving 
the DSE's These are added to the bare vertex 
ZiFjfi and the renormalization constant is fixed 
by imposing that the 7 M component is unity at a 
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renormalization point of 2 GeV. 

The non-Abelian diagram, depicted in figure Q 
dominates over the Abelian one nominally by a 
factor N%. We have checked that 1) this dom- 
inance by an order of magnitude remains after 
loop integration in the vertex for various kine- 
matics and 2) the impact of the Abelian vertex 
on the quark DSE is smaller by the same factor 
by performing a kinematic average with the ker- 
nel of the DSE as weighting function. Therefore 
to a precision of 10 % or even better one can ig- 
nore the Abelian diagram. 

2. Numerical results 

The vertex is projected into the tensor basis 
from appendix A in The loop integral is calcu- 
lated numerically in four dimensions with a stan- 
dard Gauss-Legendre grid. In one computer code 
we perform the spin sums numerically, in an al- 
ternative calculation we employ a form code to 
reduce the kernel analytically to relatively com- 
plex but tractable scalar integrals. 

If the internal qqg vertices are taken as bare 
ZifIh then the construction, that qualitatively 
has the right behaviour, is not strong enough to 
reproduce lattice data and trigger chiral symme- 
try breaking in the quark DSE. This, as com- 




Figure 1. The non-Abelian one- loop correction to 
the quark and gluon vertex is the basis of our non- 
perturbative model. To this end we replace in- 
ternal propagators by their dressed counterparts 
and enhance the internal vertices with one ghost 
dressing function each reflecting the STL 



2.5 r 

■- 



m = 115MeV; (p,) =(p/=:(p) 



X (lattice) 
(4p 2 ) X 2 (lattice) 
(-2p) X 3 (lattice) 
Xj (DSE) 
(4p 2 ) X 2 (DSE) 
(-2p) X (DSE) 




1 2 
P (GeV) 



Figure 2. Comparison of our results with lattice 
data in the particular kinematic section pi = p2- 
m q (2 GeV) is set to 115 MeV but similar results 
are obtained at the other lattice data set available 
at 60 MeV. 



merited above, is expected and can be remedied 
by enhancing the internal vertices by a ghost 
dressing factor. In figure [2] we compare the re- 
sulting qqg vertex with lattice data at the so 
called "asymmetric point" characterized by p g = 
0, pi = pi where the gluon momentum vanishes, 
Pi flows into the vertex and pi exits. This com- 
parison is successful for the leading Dirac ampli- 
tude Ai and the scalar amplitude, A3. The am- 
plitude 4p 2 A2 vanishes in our model at low mo- 
menta, whereas the lattice data (with large er- 
rors) seem to approach a constant value implying 
a divergence of A2. 

We now investigate a more interesting sec- 
tion of kinematic space, that we denote "totally 
asymmetric" point, characterized by the relations 
P2 = 2pi and p g = 3pi between the moduli of the 
momenta. This point is interesting because the 
tensor basis used is non-singular and all twelve 
different Dirac amplitudes Ai_4, Ti_g contribute 
to the vertex. We plot the four leading structures 
in figure 

The other eight Dirac amplitudes are increas- 



ingly smaller, down to two orders of magnitude 
below the plotted ones. Thus there is a rich hyrar- 
chy of Dirac amplitudes that can help in model 
building. 

3. Mass dependence and Chiral Symmetry 
Breaking 

Once we are in possession of a construction 
that successfully compares to lattice data, we 
employ it to perform a study of the mass de- 
pendence of the vertex. First note the Abelian 
one-loop diagram contains two quark propaga- 
tors, and is therefore suppressed as 1/M g in the 
heavy quark limit, whereas the non-Abelian dia- 
gram will damp as 1/M q . Therefore future cal- 
culations relating observables in the charmonium 
and bottomonium systems (see |Hj for discussion) 
will be sensitive to what class of diagrams enters 
the vertex model. 

An interesting observation is that for the range 
of quark masses considered in the lattice calcula- 
tions, the A3 Dirac amplitude has a maximum. 
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Figure 3. Leading tensor structures in the kine- 
matic section given by p\ :— p, pi — 2p, p g = 3p 
in our model construction. 



Figure 4. Mass dependence of the quotient Xi/A 
that equals 1 in an Abelian theory (as incorpo- 
rated in the Ball-Chiu construction). 



If the current quark mass is further increased, 
the intermediate quark propagator suppresses the 
vertex loop. Conversely, approaching the chiral 
limit, the quark mass function takes its minimum 
value dictated by chiral symmetry breaking alone 
and the scalar part of the loop construction also 
has a minimum (slight corrections to this result 
are expected after self-consistently solving a ver- 
tex equation). 

As for the leading Ai amplitude, we plot in fig- 
ure^ the quotient \\/A, where deviations from 
unity signal departures from the Abelian Ward- 
Takahashi identity. 

4. Outlook 

Results similar to ours have been indepen- 
dently obtained [H] in a different scheme. This 
exploits the STI for the three-gluon vertex to 
model it. Since this amounts to the resummation 
of a totally different class of diagrams we would 
need better lattice data to distinguish both mod- 
els. Our construction is of course valid (within 
approximations) for all possible kinematics and 
not just when the gluon momentum vanishes. On 
the positive side, both works concur in predicting 
p 2 \2(p 2 ) to vanish as p — -> 0. We should note that 



the existing lattice data, with large error bands, 
suggests instead this limit is finite, implying a di- 
vergence in the vertex, in the Dirac amplitude A2. 

When multiplied by appropriate powers of p 
the resulting dimensionless Dirac amplitudes of 
our vertex construction vanish as p — > except 
the leading structure Ai that takes a finite value. 

Whether a divergence can arise as a conse- 
quence of the feedback of the obtained vertex 
model on the loop construction itself (implying 
a self-consistent solution is needed) or as a conse- 
quence of the backreaction on the quark SDE is a 
topic under current scrutiny. The construction of 
a quark scattering kernel based in our model ver- 
tex and consistent with chiral symmetry is now 
straight-forward. 

5. Questions from the audience. 

• How are Euclidean space singularities in the 
two-point functions mapped to Minkowski 
space after analytical continuation? 
This is a nonperturbative problem and we 
lack a full answer at this stage. One has 
obtained some understanding of the analyt- 
ical structure of the two-point functions in 
a recent work TOI . Also Hamiltonian calcu- 
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lations in Coulomb gauge, not in Euclidean 
space, provide clear evidence for the picture 
of ghost enhancement and gluon suppres- 
sion at low momentum 

• The running a s you are employing seems 
to be somewhat too high in the middle- 
momentum range. What is the impact of 
this in the results reported? 
This study does not aim at precisely pin- 
pointing details of the propagators and ver- 
tex functions, but to obtain insight into 
their qualitative features and structure, and 
identify possible divergences that may oth- 
erwise escape lattice calculations. In this 
respect, the impact of this few-percent de- 
viation can be ignored. 




Figure 5. Exact vertex DSE equation. 



• Why is the dressing of the three-gluon ver- 
tex omitted? 

As can be observed in figure El all dressing 
of the triple gluon vertex can be absorbed 
in the quark-gluon scattering kernel. Our 
model based on the one-loop correction to 
the qqg vertex can be also viewed as an ap- 
proximation to this kernel. In this exact 
equation we are neglecting completely 
the two last terms, involving the ghost- 
quark scattering kernel and the three-gluon- 
quark scattering kernels, as their skeleton 
expansion starts at two loops. 



We thank M. Bhagwat, C. D. Roberts and P. 
Tandy for valuable discussions. The lattice data 
in figure are a courtesy of the authors of JT^. 
F. L. E. thanks the hospitality of the Tubingen 
Institute as well as a DAAD stipendium and 
Univ. Complutense travel grant. Partial support 
from grants FPA 2000-0956, BFM 2002-01003 
(MCYT, Spain), Al 279/3-4, Fi 970/2-1, and 
GRK683 (DFG, Germany). 

REFERENCES 

1. R. Alkofer and L. von Smekal, Phys. Rept. 
353 (2001) 281 |arXiv:hep-p h/0007355|. 

2. P. O. Bowman et al., |arXiv:hep -lat/0402032 



3. C. S. Fischer and R. Alkofer, Phys. Rev. D 
67 (2003) 094020 arXiv:hep-ph/0301094 . 

4. C. S. Fischer, F. Llanes-Estrada and 
R. Alkofer, |arXiv:hep-ph/ 0407294 

5. A. I. Davydychev, P. Osland and L. Saks, 
Phys. Rev. D 63 (2001) 014022 

6. J. Skullerud and A. Kizilersu, JHEP 0209 
(2002) 013 [arXiv:hep-ph/0205318| . 

7. J. I. Skullerud, P. O. Bowman, A. Kizilersu, 
D. B. Leinweber and A. G. Williams, JHEP 
0304 (2003) 047 |arXiv:hep- ph/03031 76 . 

A. Holl, A. Krass- 
and P. C. Tandy, 



10 



11 



Bhagwat, 



M. S 

nigg, C. D. Roberts 
|arXiv:nucl-th/0403012| 
M. S. Bhagwat and 
|arXiv:hep-p'h7 0407163 . 
R. Alkofer, W. Detmold, 
and P. Maris, Phys. Rev 
arXiv:hep-ph/0309077 . 
A. P. Szczepaniak and 
son, Phys. Rev. D 65 



P. C. Tandy, 



C. 
D 



S. 
70 



Fischer 
(2004) 



E. S 

(2002) 



D. 



arXiv:hep-ph/0312254 


C. Feuchter 


H. Reinhardt, 


arXiv:hep-th/0402106 



Swan- 
025012 
Zwanziger, 
and 
id. 



work in preparation. 
12. W. J. Marciano and H. Pagels, Phys. Rept. 
36 (1978) 137; E. Eichten and F. Feinberg, 
Phys. Rev. D 10 (1974) 3254. 



